Instanton analysis of hysteresis in the 3D Random-Field Ising Model 
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We study the magnetic hysteresis in the random-field Ising model in 3D. We discuss the disorder 
dependence of the coercive field H c , and obtain an analytical description of the smooth part of the 
hysteresis below and above H c , by identifying the disorder configurations (instantons) that are the 
most probable to trigger local avalanches. We estimate the critical disorder strength at which the 
hysteresis curve becomes continuous. From an instanton analysis at zero field we obtain a description 
of local two- level systems in the ferromagnetic phase. 

PACS numbers: 75.10.Nr, 75.60.Ej, 05.50.+q 
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In disordered systems whose pure counterparts un- 
dergo a first order phase transition as an external pa- 
rameter (e.g., field or pressure) is tuned, the presence 
of randomness leads to interesting hysteretic phenom- 
ena 0. Those arise from local modifications of the spin- 
odal point due to disorder, as well as from the presence of 
many metastable states. The simplest system exhibiting 
such phenomena is the random field Ising model (RFIM), 
often studied as a representative model for hysteresis in 
the presence of random fields 00000,0. 

In the RFIM, the gradual transition between states 
of negative and positive magnetization proceeds via 
avalanches triggered by the increase of the external mag- 
netic field. Previous studies 0,0] have shown that in the 
presence of strong disorder all avalanches are bounded 
and the hysteresis loop is macroscopically smooth. In 
weak disorder, however, there is a finite coercive field be- 
yond which the ferromagnetic coupling induces macro- 
scopic avalanches that span large parts of the sample 
and lead to a sharp jump in the hysteresis curve. The 
two disorder regimes are separated by a critical point 
which is believed to control the power law distribution 
of Barkhausen noise 0, EI- Recently, indications for 
such a disorder induced transitions have been reported 
in disordered magnets 0| as well as in the context of 
capillary condensation in aerogels . 

The hysteresis in the RFIM is intimately tied to the 
complexity of its free energy landscape, possessing an ex- 
ponential number of metastable minima. However, the 
precise connection between the disorder induced energy 
landscape and the nature and extent of the avalanches 
has not been established. So far, studies of the hystere- 
sis in the RFIM concentrated on lattice simulations at 
zero temperature [a , and most analytical insight was 
restricted to the Bethe lattice |2| or the extreme mean 
field limit where all spins interact with each other. 
In this Letter, we make a first step towards an analyt- 
ical theory for 3D systems. In particular, we analyze 
the features of the hysteresis loop as a function of disor- 
der strength and sample size and characterize the typical 
avalanches occurring along the hysteresis loop. This pro- 
vides a new perspective on the above-mentioned disorder 



induced transition. Further, by studying typical config- 
urations that trigger avalanches at zero field we obtain 
insight into metastability in the ferromagnetic phase. 

Before giving the details of our analysis, we summarize 
the qualitative picture we have obtained, concentrating 
on the raising branch of the saturation hysteresis loop 
(see Fig.l). In a pure system, hysteresis arises due to the 
metastability of the negatively magnetized state in ex- 
ternal fields H below the spinodal field H sp . This simple 
scenario is strongly modified by disorder: Due to rare 
regions with strong positive fields, bubbles of positive 
magnetization are created spontaneously at any small H . 
However, they are prevented from spreading by disorder 
induced pinning. The latter leads to a finite coercive field 
H c which grows with disorder strength. For fields be- 
low H c , the magnetization curve is determined by small 
avalanches, triggered by the increase of H . At H c sample 
spanning avalanches occur, while beyond H c only small 
regions with strong negative random fields may resist 
the invasion of the positively magnetized phase, finally 
collapsing under further increase of the field. With in- 
creasing disorder strength, the local avalanches become 
so dense that the emerging bubbles already percolate be- 
fore the threshold for unlimited bubble spread, H Cl is 
reached. In this situation, there is no room for avalanches 
to spread at H Cl and the hysteresis loop becomes contin- 
uous. 

In order to study the nature of avalanches in the 3D 
RFIM, we have analyzed its continuum version, described 
by the standard Landau- Ginzburg free energy 



F = Fn dx 



(yip) 2 + V{<p) -Hep- h(x)(p 



(1) 



Here <p(x) is the coarse-grained magnetization whose ten- 
dency towards ferromagnetism is described by the poten- 
tial V{ip) = ^(ip) 2 + f (v) 4 with the reduced tempera- 
ture m 2 cx t = —(1 — T/T c ). H is the uniform external 
magnetic field, and h(x.) is a Gaussian random field with 
zero mean and variance (h(x.)h(x')) = A 6(x — x'). The 
prefactor Fq depends on the microscopic details and in- 
creases with the range of ferromagnetic interactions. 
We only consider T < T c (to 2 < 0) which is neces- 
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FIG. 1: Top: Hysteresis curves for increasing disorder. At the 
coercive field H c the magnetization jumps discontinuously. In 
strong disorder (A > A c w As) the jump disappears since 
bubbles are already dense at Hd P - Bottom: The shaded areas 
indicate where local avalanches are observed in the H — A 
plane (in finite volumes Vi < V% < oo [lighter shaded] , and for 
V = oo [dark area, H < Hd P = aA]). Macroscopic avalanches 
occur at the coercive field H c (thick blue line). The horizontal 
lines correspond to the hysteresis curves in the top panel. 



sary to ensure the presence of metastable states. It is 
convenient to rescale all quantities according to natural 
units, r = x/4 *(r) = <p(x)/<po, h(r) = (£ 2 ft(x))/(2^ ), 
and H = (e 2 H)/(2tp a ), where <f a = {\m 2 \/g) 1/2 is the 
equilibrium magnetization of the pure system and I — 
(2/m 2 ) 1 / 2 ~ t~ x / 2 is the mean field correlation length. 
The reduced disorder strength is A = (A£)/(4c/?o) ~ 
(A)t~ 3 / 2 . Further, we denote by ^±{H) the minimum 
with positive/negative magnetization of V(^) — H^>, 
where V(W) = -\^ 2 + j* 4 . 

In order to study the free energy landscape of this 
model, we restrict ourselves to a mean field analysis and 
concentrate on local extrema of Q , satisfying 



(Stf(r) 2 { ' 



V"(*(r)) - H - h(r) = 0. (2) 



In the spirit of the rate independent hysteresis approxi- 
mation [13 we focus on spontaneously occurring events 
and neglect thermal activation. This is justified suffi- 
ciently far below T c where typical free energy barriers 
between local minima are larger than T |14| . 

We may restrict the analysis of the hysteresis loop to 
the raising branch where the magnetization is everywhere 



close to ^/-(H) at sufficiently negative H. As H is in- 
creased, the magnetization either remains in its local free 
energy minimum, responding paramagnetically to the 
field change, or it becomes locally unstable, and a bub- 
ble with higher magnetization is created spontaneously. 
This requires a sufficiently strong local excursion h(r) of 
the random field. For weak disorder such configurations 
are rare, their density scaling as p oc exp[— S], where 
S = J dr[/i(r)] 2 /2A. Hence, we may think of them as 
dilute islands each of which can be analyzed separately. 
Further, we may assume that at some distance from a 
rare excursion, the random fields are weak such that the 
magnetization tends to 

The fate of a spontaneously emerging bubble depends 
on its size. We will see that typical bubbles are small and 
do not spread much since they are held back by a surface 
tension barrier. However, this does not apply to larger 
(and hence exponentially rarer) bubbles. Instead, their 
dynamics is governed by the competition between the ex- 
ternal field exerting a driving force 2H on the surface of 
the bubble, and the random field disorder, which tends 
to pin it. Collective pinning theory predicts the critical 
field for domain wall depinning to scale as H^p = a A |l5j . 
From a high velocity expansion [lfi| we obtained the es- 
timate a « 0.055, which we confirmed by numerical sim- 
ulations [l7j . For H < H^p, all avalanches are bounded, 
and the hysteresis loop is dominated by typical, i.e., most 
abundant, small avalanches. Beyond -ffdp, the first un- 
stable bubble which is not constrained by surface tension 
induces a sample spanning avalanche and a macroscopic 
jump in magnetization, however low the probability per 
unit volume for such an event may be. In the thermody- 
namic limit, the coercive field H c thus coincides with the 
depinning field H^p- We will discuss mesoscopic fluctu- 
ations of H c in weakly disordered small samples further 
below. The approach to complete saturation beyond H c 
is governed by small bubbles of negative magnetization 
that spontaneously implode upon an increase of field. 

Let us now analyze in more detail the properties of typ- 
ical avalanches as a function of disorder and field, both 
for H < H c and H > H c . An avalanche is triggered spon- 
taneously when a local free energy minimum (with mag- 
netization ^(r)) becomes unstable upon increase of H. 
This occurs when \l/(r) merges with an adjacent saddle 
point, 4 , sp (r), the magnetization difference S^> = \& sp — \& 
turning into a soft mode of the Hessian, 



-1/2V 2 + V r "(*))5* = 0. 



(3) 



Let us denote by M. the set of all states such that 
there is a soft mode 5^> satisfying Eq. 10). These states 
are "marginal" in the sense that they are unstable with 
respect to the smallest perturbations. 

Typical avalanches are triggered by the most proba- 
ble configurations ( "instantons" ) of random fields, h — 
h>iast( r \H'), such that the current magnetization state, 
^(r), becomes marginal (\? e M.) when the external 
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field reaches H' . These instantons thus correspond to 
non-trivial local minima of the action 



5 =2A 



2 v ' 



i? 



(4) 



within the set of marginal states W e A4. Here, we used 
Eq. 101 to express the random field in terms of VF Assum- 
ing spherically symmetric instantons, the Euler-Lagrange 
equation takes the form 



-AW 



2 

A (d+2) 



+ V"(*) 
+ V"(¥) 



A 1 J * T// 



7'(¥)+JT EE (5) 



V'(tf) + £n = Cn* 



where A^ d ' = <9 2 + ^-^-d r denotes the radial part of the 
Laplacian in d dimensions, C is a Lagrange multiplier, 
and nq, is the local normal to the manifold M.. 

Performing an extensive search for local minima of S 
(analytically by solving J^J and numerically by varia- 
tional minimization of Q) we found two types of optimal 
random field configurations. For H < Hi = 0.0763, the 
optimal random field is such that the state \& corresponds 
to a threefold degenerate free energy minimum. However, 
slight deviations 6h(r) from the instanton configuration 
/linst lead to a local free energy landscape where the state 
with lowest magnetization is marginally unstable with 
respect to a nearby minimum with locally higher mag- 
netization (5M cc (j Sh 2 ) 1 / 4 ). An essentially identical 
analysis for H c < H < +oo, but with inverse boundary 
conditions — > yields similar instantons corre- 

sponding to the dominant bubbles which implode upon 
increase of the external field beyond H c (lif. 

A second type of instantons describes the typical 
avalanches in the interval H\ < H < H c . These instan- 
tons can be found analytically by solving the "dimen- 
sionally reduced" equations of motion, -1/2 A' d " 2 '$ + 
V'(^t) — H = 0, which indeed yield a solutions of 
Eq. (jnj with C = 0. The corresponding random field 
is hi nst (r\H) = —(d r ty)/r. In d = 3 dimensions, this 
equation is readily integrated with boundary conditions 
<T(r = 0) = and f(Y -> oo) -> ®-(H). The action 
Eq. Q associated to these field configurations decreases 
monotonically from S(Hi) = 5.7/A to S(H sp ) = 0. 

An analysis of the local free energy landscape associ- 
ated to this type of instantons shows that for H i < H < 
H 2 = 0.173 the typical bubbles are constrained by a sur- 
face tension barrier, the change of magnetization remain- 
ing local (see Fig. |2Ja)). This is illustrated in the lower 
panel of Fig. [21 where we plot the instanton hi nst (r\H) 
together with all solutions of Eq. © for H = 0.13. Ini- 
tially, the system is in the marginal state (I) with the 
lowest magnetization. Upon slight increase of H, it be- 
comes unstable and the system runs to the nearby lower 
lying minimum (II) which represents a localized bubble 
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FIG. 2: Top: Schematic sketch of the free energy landscapes 
associated with typical avalanches for various regimes of the 
external field H. (a) H < H2 (cf. bottom panel): The state 
(I) becomes unstable and runs to the metastable (II) with a 
local bubble of higher magnetization. It is separated from 
the global minimum (IV) by a surface tension barrier (III); 
(b) H = H2: state (II) and (III) merge and the surface barrier 
vanishes; (c) H > H2 '■ Rare bubbles in negative random fields 
collapse under further increase of H . Bottom: Solutions (I- 
IV) of Eq. obtained numerically for H = 0.13 < H2 and 
h = hi nst (r\H) [thick line]. 



of positive magnetization. Further spread of this bubble 
is prevented by the surface tension of the domain wall 
that separates the positively and negatively magnetized 
regions. Indeed, a further saddle point (III) separates 
the local minimum from the global minimum (IV). Note 
that as long as H < -ffdp, domain wall pinning by typ- 
ical random fields will split the state (IV) into a set of 
metastable states which differ on larger scales which are 
not shown in Fig. [21 

Beyond H2, the confining surface tension barrier 
vanishes, and the avalanches spread further out (cf. 
Fig. [2{b,c)), leading to a marked jump in the slope of 
the hysteresis curve. However, this effect is only observ- 
able in a small window of disorder (see Fig.^J. Indeed, 
for A < H2/ct, macroscopic avalanches set in before H2, 
while for A > H2/a the density of bubbles becomes 
rapidly high, preventing large avalanches as discussed be- 
low. 

At a critical disorder A c the typical bubbles start to 
percolate at the depinning field f/d P (see Fig. 0). This 
marks the crossover towards a regime with a continu- 
ous hysteresis loop. Indeed, for A > A c bubbles are 
dense already below the depinning threshold Hd P , such 
that avalanches at H^p run into already positively mag- 
netized regions and thus die out. We can estimate A c 
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FIG. 3: Schematic view of the temperature evolution of the 
local free energy landscape in a random field hi ns t(r\H — 0) + 
Sh. Upon cooling, the high temperature state (H) becomes 
unfavorable with respect to an emerging low energy state with 
more uniform magnetization (L), but remains metastable in 
a temperature range AT due to a free energy barrier AT. 



from the requirement that the density of typical bub- 
bles is p(H dp [A c }) = O(l), or S(H dp [A c ]) ~ 1. A more 
careful calculation including Gaussian prefactors yields 
A c ~ 3.2, which translates into A c ~ 9 (\m 2 \) 3 / 2 /g ~ 
const x r 3 / 2 / 1 g for the transition line in the T — A plane. 
It is interesting to note that for A « A c the depinning 
field iidp = aA is very close to the threshold field H2 for 
typical avalanches to overcome surface tension barriers. 

Above we used the identity of the coercive field and the 
depinning threshold. This applies to the infinite volume 
limit, as a result of which one finds H C [A — > 0] — > (see, 
e.g., 0), instead of H C [A — > 0] — > H sp , which holds for 
finite size systems (see bottom of Fig. 0, H sp = 2/3v3 w 
0.385 being the spinodal field of the pure system. For 
a finite system of volume V, the only type of bubbles 
occurring in the hysteresis loop are those with sufficiently 
large density, Vp(H,A) > 1. This relation allows us to 
estimate the field H* at which the first (typical) bubbles 
occur. If the disorder is so weak that H* > H 2 , there 
will be a single sample spanning avalanche around the 
coercive field H* ~ H sp — 0.5Alog[V/£ 3 ]. For stronger 
disorder (H* < H2) local bubbles occur at H* while only 
at higher field a macroscopic avalanche is triggered by the 
first bubble that is not constrained by surface tension. 

The analysis of instantons also provides information on 
the emergence and bifurcation of metastable states upon 
temperature variations. We have studied the tempera- 
ture evolution of the free energy landscape in the presence 
of local random field configurations close to the instan- 
tons corresponding to H = and T — T , h = h inst (H = 
0; To) + Sh. These describe the rare regions which locally 
admit two metastable configurations in the ferromagnetic 
phase, cf. Fig. [3] At high temperatures there is only one 
state (H) in which 'F(r) adjusts to the random field. At 
lower temperature, a second state (T) with more uniform 
magnetization emerges and eventually dominates. This 



confirms a scenario for the emergence of metastability 
conjectured by Villain in the context of pinned domain 
walls 0|- Over a range of order AT ~ 5 3 ^ 2 (where 
6 = (J Sh 2 ) 1 / 2 ) the two states coexist, being separated 
by a free energy barrier of order AT ~ 5 2 . The spatial 
density of such two-level systems with a given free energy 
barrier scales as p(AF) cx exp [— So — 0(AF 1 ^ 2 )\ with 
So = S[h inst (r\H = 0)] ~ 8.2/A. These two-level sys- 
tems are very similar to the metastable bubbles observed 
numerically above the ferromagnetic transition |2flj . We 
expect a deeper analysis of the bifurcation of states along 
the lines of this Letter to yield new insight into the nature 
of the phase transition in the RFIM. 
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